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Abstract: 

_ ^ 

■  A  seismic  refraction  experiment  was  conducted,  as  part  of  the 
November-December,  1973  field  tests  off  Daytona  Beach,  Florida.  In 
the  seismic  experiment  reported  here,  two  pound  charges  were  detonated 
one  meter  off  the  bottom  and  received  on  a  vertical  hydrophone  array. 
The  experiment  was  conducted  along  two  legs;  a  constant  depth  leg  and 
a  leg  along  a  sloping  bottom.  Along  each  leg,  charges  were  detonated 
at  ranges  from  about  1  km  to  10  km.  The  recorded  signals  along  the 
constant  depth  leg  have  been  examined  and  four  distinct  refracted 
arrivals  have  been  identified.  The  travel  time  versus  distance  curves 
for  these  different  arrivals  have  been  analyzed  and  layer  thicknesses 
and  propagation  velocities  for  the  different  sedimentary  layers  have 
been  obtained.  The  data  presented  here  gives  a  description  of  the 
bottom  properties  that  are  important  in  modeling  the  sound  propagation 
characteristics  for  this  area. 
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Introduction 

During  the  field  trip  conducted  in  November  and  December,  1973  in 
conjunction  with  NRL,  a  seismic  refraction  experiment  was  conducted.  In 
this  experiment,  two  pound  charges  were  deployed  from  the  ship  moving  at 
six  knots,  and  detonated  one  meter  from  the  bottom.  The  signals  were 
received  on  a  vertical  hydrophone  array  and  transmitted  to  the  ship  where 
the  signals  were  recorded  on  magnetic  tape.  The  amplifier  gains  were 
adjusted  so  that  the  refracted  signal  arrivals  in  the  bottom  (ground  waves 
could  be  detected.  The  shots  were  made  along  two  legs;  a  constant  water 
depth  leg,  and  a  leg  along  a  sloping  bottom,  with  the  propagation  from 
shallower  to  deeper  water.  On  each  leg,  the  shots  were  made  from  ranges 
of  about  1  km  to  10  km.  In  this  report,  the  results  of  the  constant  depth 
leg  experiment  will  be  examined. 

A  hydrophone  (instant  shot  hydrophone)  was  trailed  about  200  feet 
behind  the  ship  and  the  signal  received  from  the  explosion  on  this  hydro¬ 
phone  provided  a  time  base  for  determining  the  travel  time  for  the  various 
signal  arrivals  on  the  vertical  receiving  array.  The  results  to  be  pre¬ 
sented  here  were  received  on  the  bottom  hydrophone,  which  was  about  three 
meters  from  the  bottom. 

Analysis 

Since  the  instant  shot  hydrophone  was  not  located  at  the  point 
where  the  charge  was  detonated,  we  must  make  a  correction  to  the  travel 
time  determined  by  the  time  difference  between  the  instant  shot  time  and 
the  time  of  arrival  of  the  various  signals.  The  geometry  is  shown  in 
Figure  1. 

The  ship  was  moving  at  six  knots  and  the  charges  had  a  65  second 
fuse.  Thus  d,  the  distance  from  the  stem  of  the  ship  to  the  point  where 
the  shot  is  detonated,  is  given  by 

d  -  (6  knots)  (65  sec.)  =  658.7  feet. 


(1) 
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FIGURE  1 


The  water  depth  was  about  42  meters  and  the  charges  were  detonated  one 
meter  from  the  bottom.  The  instant  shot  hydrophone  was  streaming  200 
feet  behind  the  ship,  and  near  the  surface.  We  will  assume  that  the 
instant  shot  hydrophone  was  120  feet  above  the  charge  depth.  Thus, 

Ar  =  [(658.7  -  200)2  +  (120)2]1/2 

=  474  feet  (2) 

=  144.5  meters. 


Hence,  since  the  measured  sound  velocity  in  the  water  was  1533  m/sec., 
the  time  correction  is  given  by 


At 


144.5  m 
1533  m/sec. 


=  0.094  sec. 


(3) 


In  Figures  2  through  11,  the  received  signals  for  the  ten  shots  on 
the  constant  depth  leg  for  the  bottom  hydrophone  are  shown.  These  repro¬ 
ductions  do  not  always  show  the  different  signal  arrivals  as  clearly  as 
the  raw  data  recordings  nor  is  some  of  the  supplementary  data  shown  here. 
Notice  also  that  the  time  scale  is  not  the  same  on  all  figures.  However, 
if  we  examine  these  figures,  we  see  that  as  many  as  four  signal  arrivals, 
in  addition  to  the  water  wave  signal,  can  be  seen.  The  travel  times  for 
each  of  these  various  signals  are  given  in  Table  1,  where  the  correction 
given  in  (3)  has  been  made.  The  ranges  for  the  various  shots  were  deter¬ 
mined  from  the  time  of  arrival  of  the  water  wave  signal  and  the  known  sound 
velocity  in  the  water.  The  numbers  in  parentheses  in  Table  1  will  be 
explained  later. 


Travel  Time ,  x ,  seconds 
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The  travel  times  versus  ranges  have  been  plotted  in  Figure  12  and 
we  see  that  the  arrivals  lie  nearly  on  straight  lines.  If  we  find  the 
best  linear  fit  to  minimize  the  mean  square  error,  we  obtain 


E  t  R.  -nr  R 
1  1  1 

T  -  T  = 


J  (V2  -  n(R)2 


(R  -  R) 


(4) 


where  the  overbar  denotes  average  value.  Using  (4),  we  find  that  the  best 
linear  fits  are  given  by 


Line 

A: 

T 

=  0.0005804R  + 

0.00122 

(5) 

Line 

B: 

T 

=  0.0005308R  + 

0.01801 

(6) 

Line 

C: 

T 

=  O'.  0004403R  + 

0.1339 

(7) 

Line 

D: 

T 

=  0.0002617R  + 

0.5417 

(8) 

Substituting  the  ranges  given  in  Table  1  into  these  equations,  we 
find  the  values  of  travel  time  shown  in  parentheses  in  Table  1.  Notice 
that  the  agreement  is  excellent  and  that  the  points  all  lie  very  nearly 
on  the  computed  lines.  Note  also  thac  the  arrival  at  0.633  seconds  at 
a  range  of  1,126  meters  should  probably  have  been  assigned  to  line  C 
rather  than  line  B . 


We  now  want  to  relate  the  observed  travel  time  versus  range  data 
to  the  bottom  parameters.  We  will  assume  a  layered  bottom  with  horizontal 
interfaces. 


Source 


c^  ■  1533  m/sec. 


Receiver 
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Figure  13  shows  the  geometry  for  arrival  A.  The  sound  hits  the 
bottom  at  the  critical  angle,  is  refracted  along  the  interface  and  travels 
with  the  velocity  C2*  It  leaves  the  interface  at  the  critical  angle  and 
travels  to  the  receiver.  We  have  that 
C1 

sin0,  =  — 

1  c2 


11  - 


c2  ~  C1 


=  r^sine^  = 


/  2  2' 
(2  "  C1 


/C  2  -  C1 


lc2  ~  C1 


D1  =  R  "  X1 


Then,  the  travel  time,  x,  is  given  by 


rl  D1  r2  °1 

C1  c2  C1  c2 


rl  +  r2 


rl  +  r2 


X1  +  x2 


FIGURE  14 
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We  have 


sin0o  =  — 

3  c3 


r3  = 


T  = 


c3H3 

JC  3  - 

- \ 

2 

C2 

C2H3 

\|c  3  " 

2 

c2 

D 

JL  + 

2r3 

c3 

c2 

R  . 

2h3/ 

c3 

we  have 

1 

C2C3 


“3  0.0005308 


2“  Jl]-  4' 

and  — JJLi - £ —  =  0.01801. 


C2C3 


Thus , 


=  1884  m/sec.  and 


=  38.4  m. 


We  will  now  examine  arrival  C  as  shown  in  Figure  15 . 
Source 


T~X 

h, 


\ 


(11) 

(12) 


Receiver 
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We  have 

sin0/ 


sin0r 


sin0c  =  — 

5  c, 


Using  elementary  geometry,  as  above,  we  find 


r  ,  2hH-4'  , 

T  =  -  +  - -  + 

C4  C4C2 


2HWc4-g3 
C4C3 


But,  from  (7)  we  have 

x  =  0.0004403  R  +  0.1339. 

Thus , 

=  2271  m/s ec. 


0.0004403 


and 


2H3  v  4  "  f2_  +  “*‘4 


2H . 


/  2  2' 

Vc4  c3 


C4C2 


C4C3 


=  0.1339. 


Using  the  values  found  for  C2,  c^,  and  in  (9),  (1),  (14)  and 
we  obtain, 


H.  =  176.9  m. 

4 

Finally,  we  will  consider  arrival  D  as  shown  in  Figure  16. 
We  have  that 

sin0,  sin9_  sin0„  c. 

- i  =  - 1  =  - §_  >  gln0  =  . 

c2  c3  c4  8  c5 


1975 

R-1213 


(13) 


(14) 


(12), 


(15) 
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Source  Receiver 


FIGURE  16 


By  the  same  procedure  as  before,  we  find  for  the  travel  time, 


t  =  —  + 

C5 


Z»J4  -  c* 


■\ 


C2C5 


But  from  (8) ,  we  have 

x  =  0.0002617  R  +  0.5417. 

Thus , 

=  3820  m/sec. 


C5  “ 


0.0002617 


2H4  V  C5  ~  C3 


C3C5 


2H5  v  5  " 


C4C5 


(16) 


and 


2H3  fCl  -  1  +  2RJCl  ~  C3 


2H5  /C5  ~  C4 
C4C5 


C2C5 


C3C5 


=  0.5417. 


Using  the  values  for  c2>  c3>  c^,  H3  and  found  above,  we  obtain 

(17) 


=  478.1  m 
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Discussion  of  Results 

Based  on  the  analysis  given  above,  we  postulate  that  the  sedimentary 
structure  in  the  area  where  the  seismic  refraction  experiment  was  made  is 
that  shown  in  Figure  17.  In  reference  [1],  the  results  of  a  large  number 
of  seismic  refraction  profiles  are  given  for  the  same  general  area  as  the 
profile  analyzed  here.  In  Table  2  of  that  report,  layer  velocities  and 
layer  thicknesses  are  given  for  profiles  made  in  continental  shelf  loca¬ 
tions.  For  the  seismic  profiles  that  were  reported  in  that  paper,  much 
longer  ranges  and  much  larger  charges  were  used  than  for  the  profile 
reported  here.  Consequently,  information  was  obtained  for  much  deeper 
layers  than  was  obtained  from  our  data.  On  the  other  hand,  since  the 
profile  was  obtained  by  using  short  range,  closely  spaced  shots  in  the 
experiment  reported  here,  the  resulting  data  on  the  uppermost  sedi¬ 
mentary  layers  should  be  more  reliable  than  that  reported  in  reference  [1]. 

If  we  examine  Table  2  of  reference  [1],  we  see  that  a  low  velocity 
layer,  denoted  by  Layer  A,  was  observed  on  nearly  all  the  profiles.  The 
velocity  associated  with  this  layer  varied  from  1660  m/sec.  to  1900  m/sec., 
with  layer  thicknesses  ranging  from  50  m  to  580  m.  The  velocity  found 
for  the  uppermost  layer  in  the  test  reported  here,  1723  m/sec.,  is  consis¬ 
tent  with  these  results.  However,  the  thickness  we  found,  38  m,  is 
slightly  less  than  any  reported  in  reference  [1].  The  second  layer  shown 
in  reference  [1],  Layer  B,  was  observed  on  only  a  few  profiles,  and  in 
only  one  case  was  both  Layer  A  and  Layer  B  present.  The  velocities  reported 
for  this  layer  ranged  from  1920  m/sec.  to  2050  m/sec.,  with  thicknesses 
varying  from  60  m  to  630  m.  Thus,  the  velocity  observed  here,  1884  m/sec., 
is  slightly  less  than  that  found  in  reference  [1].  The  layer  thickness 
observed  for  this  layer  is  consistent  with  the  earlier  results.  The 
next  layer,  C,  reported  in  the  earlier  work  was  the  most  consistent  layer 
since  it  was  observed  on  every  profile  except  one.  The  signal  arrival 
associated  with  this  layer  was  the  most  prominent  feature  observed  in  our 
received  signals,  i.e.  the  signal  was  very  strong  and  very  distinct.  The 
velocities  found  for  Layer  C  in  reference  [1]  varied  from  2250  m/sec.  to 
2880  m/sec.,  with  thicknesses  varying  between  190  ra  and  2.45  km.  Thus, 
the  results  we  obtained,  2271  m/sec.  and  478  m  are  consistent  with  the 
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reported  results.  The  next  layer  reported  in  reference  [1]  had  velocities 
between  3160  m/sec.  and  3880  m/sec.  The  velocity  obtained  here  for  this 
layer  was  3820,  which  again  is  consistent  with  the  earlier  results.  In 
reference  [1],  results  are  also  given  for  a  number  of  deep,  higher 
velocity  layers,  which  were  not  observed  in  our  experiment  due  to  the 
short  ranges  that  were  used. 

Recommenda  t ions 

Since  the  sediments  deeper  than  several  water  depths  are  not 
expected  to  have  a  significant  effect  on  sound  transmission,  except  for 
extremely  low  frequencies,  it  is  recommended  that  the  deeper  layers  be 
neglected  in  modeling  the  bottom  for  this  location.  Based  on  the  results 
of  this  experiment,  the  bottom  in  this  area  should  be  modeled  by  a  layer 
of  velocity  1723  m/sec.  with  a  thickness  of  38  m  on  top  of  a  half-space 
with  velocity  1884  m/sec. 
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